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I.  INTRODUCTION 


This  report  documents  the  idea  or  concept  of  operating  an  existing  Silicon-on-Insulator 
(SOI)  Micro-ElectroMechanical  Systems  (MEMS)  gyroscopic  sensor  previously  developed  in 
the  early-2000s  under  MEMS-based  Angular  Rate  Sensor  (MBARS)  and  MicroControlled  Array 
Sensors  (mCAS)  [1,2,3].  This  report  serves  as  documented  evidence  for  future  use  of  this 
open-sourced  idea  of  the  measurement  of  two-axis  of  linear  acceleration  from  a  previous  sensor 
originally  designed  to  operate  as  a  single-axis  rotation  sensor. 

II.  GYROSCOPE  CHIP  AS  A  TWO-AXIS  ACCELEROMETER 

This  section  presents  the  development  of  the  equations  that  enable  the  gyro  system  to 
function  as  a  two-axis  (X-axis  and  Y-axis)  accelerometer  while  functioning  as  a  Z-axis  rate  gyro. 
When  the  mass  of  the  gyro  is  subjected  to  a  force  due  to  vibration,  the  mass  is  displaced  in  the 
direction  of  the  force.  This  displacement  is  in  addition  to  and  independent  of  any  gyro 
oscillations  involved  in  rotation  rate  measurements,  but  motions  due  to  vibrations  can  affect  the 
gyro’s  oscillation  pattern.  If  the  frequency  of  the  vibration  is  much  less  (a  factor  of  10  or  more) 
than  the  resonant  frequency,  these  displacements  can  be  accurately  detected  by  the  gyro’s 
electronics/software.  As  the  vibration  frequency  approaches  the  excitation  frequency  ( co ),  the 
resultant  motion  appears  to  be  at  the  beat  frequency  between  the  vibration  frequency  and  co.  The 
physical  effect  detected  depends  on  the  axis  (excitation  or  sense)  in  which  the  displacement 
occurs. 

If  the  displacement  is  in  the  direction  of  the  sense  axis,  then  the  vibrations  can  be  detected 
as  the  average  displacement  of  the  mass  for  each  co  period.  The  sense  circuit  is  designed  to 
measure  the  amplitude  of  any  oscillation  produced  by  the  rotation-induced  Coriolis  force  at  co. 
The  electronics  accomplish  this  by  generating  sense  signals  at  256  times  co  (approximately  1 
megahertz)  and  measuring  the  mass’  position  256  times  each  excitation  period.  The  software 
can  then  demodulate  this  signal  using  a  method  described  in  the  appendix,  as  originally 
published  in  a  2004  contract  report  [4],  An  average  non-zero  mass  displacement  would  manifest 
itself  as  the  constant  non-zero  value  in  the  Inphase  and  Quadrature  (IQ)  demodulation  of  the 
sense  signal.  That  is,  to  measure  the  displacement  due  to  vibration,  sum  all  the  partial  signal 
sums  computed  in  the  IQ  co  demodulation.  This  calculation  is  basically  a  by-product  of  the  sense 
demodulation  that  is  already  being  computed  and  adds  little  extra  computation  time. 

If  the  displacement  is  in  the  direction  of  the  excitation  axis,  then  the  vibrations  can  be 
detected  as  the  average  difference  between  the  two  excitation  diametrically  opposite  capacitors. 
Any  mismatch  between  these  capacitors  produces  a  signal  at  co,  as  opposed  to  the  excitation 
amplitude  signal  produced  at  2  co.  Other  signal  mismatches  (that  is,  mismatch  between  bias 
voltages  and  excitation  voltage  amplitudes)  also  produce  signals  at  co;  however,  these  are 
constant  with  time  but  do  vary  with  temperature.  Therefore,  to  detect  a  vibration  in  the 
excitation  axis  direction,  the  software  needs  to  demodulate  the  excitation  signal  at  co.  The 
software  is  already  demodulating  the  signal  at  2  co  to  track  the  resonant  frequency.  To 
demodulate  the  excitation  signal  at  co,  partial  signal  sums  need  to  be  summed  in  a  different  order, 
as  explained  in  the  appendix.  This  requires  very  little  additional  computation  time. 
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Thus,  both  the  X-axis  and  Y-axis  accelerations  can  be  measured  at  the  same  time  while  the 
chip  is  functioning  as  a  gyro  with  very  little  additional  processor  computational  time  and  no 
additional  electronic  circuitry. 

III.  SENSE  AXIS  SATURATION  AND  HIGH  ROTATION  RATE  MEASUREMENTS 

This  section  presents  the  development  of  the  equations  that  describe  the  sense  axis 
saturation  phenomenon  and  the  excitation  axis  amplitude  attenuation  phenomenon.  These 
equations  also  suggest  a  method  to  measure  the  MEMS  gyro  rotation  rate  even  after  the  sense 
channel  output  has  saturated.  In  testing  gyros,  it  has  been  noted  that  as  the  rotation  rate  increases 
(1,000  to  2,000  degrees  per  second),  the  sense  channel  output  increases  in  a  non-linear  manner 
and  usually  saturates  somewhere  between  1,500  and  2,000  degrees  per  second.  In  the  past,  this 
effect  has  been  attributed  to  the  electronics  saturating  and  spring  softening,  but  this  section 
develops  equations  that  explain  the  phenomenon  and  suggest  that  the  saturation  is  due  to  cross- 
coupled  Coriolis  forces  in  the  gyro.  The  equations  suggest  a  way  to  measure  the  rotation  rate 
beyond  the  rate  at  which  the  sense  signal  stops  increasing.  They  also  suggest  methods  to 
linearize  the  sense  and  excitation  outputs  over  limited  ranges. 

In  essence,  as  the  rotation  rate  increases,  the  sense  signal  increases  but  at  an  ever- 
decreasing  rate  and  eventually  saturates.  Also,  the  excitation  signal  decreases  but  at  an  ever- 
increasing  rate  even  past  the  rotation  rate  at  which  the  sense  signal  saturates.  The  following  is  a 
simple  example  of  an  algorithm  that  uses  both  the  sense  and  excitation  signal  to  calculate 
rotation  rate: 

•  Detennine  the  break  rotation  rate  (fie),  that  is,  where  {2mCLBcoKco)  =1,  for  the  gyro 
chip. 

•  Use  this  to  linearize  the  sense  signal  up  to  a  predetermined  rotation  rate  (f2i  =  fhs/2). 

•  Use  the  sense  output  to  compute  the  rotation  rate  between  0  and  fii. 

•  Use  this  to  linearize  the  excitation  signal  from  f2i  on.  Use  the  excitation  output  to 
compute  the  rotation  rate  above  Tf. 

The  following  subsections  develop  the  equation  for  cross-coupled  Coriolis  forces  in  the 
MEMS  gyro. 

A.  Background 

There  are  two  paradigms  in  which  vibration  gyroscopes  operate — the  constant 
excitation  amplitude  and  excitation  force  paradigms.  Both  paradigms  assume  that  a  mechanical 
mass/spring  system,  free  to  oscillate  along  two-axis  (an  excitation  axis  (X-axis)  and  a  sense  axis 
(Y-axis)),  is  subject  to  a  rotation  around  the  third  axis  (Z-axis).  A  periodic  electrostatic  force  is 
applied  along  the  excitation  axis.  Rotation  produces  an  oscillation  along  the  sense  axis,  the 
amplitude  of  which  is  proportional  to  the  rate  of  rotation. 

Equations  1  and  2  are  relevant  equations  of  motion  for  this  system,  where  x  is  the 
excitation  axis  and  y  is  the  sense  axis: 
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(1) 


x  =  - co\x - —  x - —  +  xQ2  +  yti  +  2Qv  -  a  x  cos  6-  aY  sin  0 

Q  m 

and 

y  =  -a>y y  -  ~y  ~  2Qx  +  vQ2  -  xQ  +  ax  sin  0  -  aY  cos  6  .  (2) 

Equations  1  and  2  have  units  of  acceleration  (in  meters  per  second  squared)  but  can  be 
converted  to  force  by  multiplying  by  the  mass  (m)  (in  kilograms)  of  the  proof  mass. 

For  this  discussion,  the  previous  equations  will  be  simplified.  This  will  emphasize  the 
results  being  derived.  If  these  assumptions  are  not  made,  the  resultant  equations  will  be  almost 
the  same  as  the  ones  derived  but  more  algebraically  complex.  Assume  there  are  no  external 
accelerations,  that  is,  the  m  is  not  subject  to  a  gravitational  acceleration  or  external  vibrations, 
then  ay  =  ax  =  0.  Also,  assume  that  both  axes  have  the  same  resonant  frequency,  that  is, 
co  x  =  co  y  =  co,  the  resonant  frequency  of  the  mechanical  system.  Then,  these  equations  can  be 
rewritten  as 


+  2Qv  -  =  x  +  x  —  +  x(a>°  -  Q2 )  ^ 

m  Q 

and 

-2Qx-xQ  =  y  +  y^  +  yip)2  -Q2).  (4) 

The  right-hand  side  of  each  equation  is  the  equation  of  motion  for  a  mass/spring 
mechanical  system.  The  left-hand  side  of  each  equation  shows  the  forcing  functions  (as 
accelerations)  and  the  cross  coupling  between  the  two-axis. 

Further  simplify  Equations  3  and  4  by  assuming  that  the  system  is  interrogated  only  at 
constant  rotation,  that  is,  Q  —  0 .  Also,  realize  that  co  is  greater  than  3,000  hertz,  while  Q  is 
approximately  8.3  hertz  at  a  maximum  rotation  rate  of  3,000  degrees  per  second.  Note  that  (co  - 
Q2)  ~  co2.  Thus,  Equations  3  and  4  can  be  rewritten  as 


2Qv 


Fy  ...  co  2 

- . - =  X  +  x — I-  xco 

m  Q 


and 


(5) 


.  UJ  2 

-2Llx  =  y  +  y—  +  yco  . 


(6) 


Note  that  there  is  a  set  of  cross-coupled,  second  order,  linear,  differential  equations. 
Assume  that  Fx  is  a  sinusoid  at  the  resonant  frequency,  that  is,  Fx  =  Fry\n(co  t). 
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Looking  at  the  previous  equations,  note  that  when  there  is  no  Z-axis  rotation  (that  is, 

=  0),  the  equations  are  not  coupled,  there  is  no  Coriolis  acceleration,  and  the  X-axis  oscillation 
is  only  caused  by  the  electrostatic  forcing  function.  Depending  on  initial  conditions,  the 
amplitude  of  the  X-axis  oscillation  asymptotically  grows  until  the  forcing  function  energy  is 
equal  to  the  energy  being  dissipated  by  the  dissipation  constant  (heat  loss),  and  then  the 
amplitude  remains  constant. 

A  rotation  along  the  Z-axis  produces  a  Coriolis  acceleration  along  the  Y-axis  (left-hand 
side  of  Equation  6)  which  sets  up  a  sinusoidal  oscillation,  the  amplitude  of  which  can  be 
measured  and  is  proportional  to  the  rotation  rate.  The  Y-axis  oscillation  in  the  presence  of  the 
Z-axis  rotation  produces  an  oscillation  along  the  X-axis,  as  seen  in  Equation  5.  The  X-axis 
acceleration  opposes  the  electrostatic  forcing  acceleration.  If  uncompensated  for,  the  X-axis 
Coriolis  acceleration  will  reduce  the  amplitude  of  the  X-axis  oscillation. 

B.  TheCEAP 

In  the  Constant  Excitation  Amplitude  Paradigm  (CEAP),  as  the  gyro  is  rotated,  the 
amplitude  of  the  X-axis  oscillation  is  maintained  by  increasing  the  electrostatic  force’s  amplitude 
to  counteract  the  X-axis  Coriolis  acceleration.  Thus,  the  electrostatic  forcing  function  can  be 
divided  into  two  terms: 


Fx=FZ+F£,  (7) 

where,  F£  is  the  constant  electrostatic  force  needed  to  maintain  a  constant  X-axis  oscillation, 
and  F£  is  the  electrostatic  force  needed  to  just  balance  the  X-axis  Coriolis  force,  that  is: 

Fc  =  2mD.(t)y(t )  .  (g) 

Since  the  amplitude  of  the  X-axis  oscillation  is  maintained  at  a  constant  level,  the 
Y-axis  Coriolis  force  equation  can  be  written  in  a  straightforward  way  as 

Fc  =  -2mQ(t)x(t) ,  (9) 

where,  m  is  the  gyro’s  mass,  Q(t)  is  the  Z-axis  rotation  rate  (in  radians),  and  x(t)  is  the  velocity 
of  m  along  the  X-axis  (see  Equation  6).  Since  the  amplitude  of  x(t)  is  being  maintained  at  a 
constant  value  for  Coriolis  measurements,  there  is  no  need  to  be  concerned  about  the  Y-axis 
oscillations  being  fed  back  to  the  X-axis  because  it  does  not  have  any  impact  oni(t)  . 

Therefore,  in  CEAP,  the  amplitude  of  the  excitation  oscillation  must  be  measured.  This 
amplitude  measurement  is  used  in  a  closed-loop  feedback  control  system  to  continuously  change 
the  gain  of  the  forcing  function  amplifier  as  the  rotation  rate  changes  in  order  to  maintain  a 
constant  X-axis  oscillation  amplitude.  In  an  ideal  system,  the  amplitude  could  be  perfectly 
maintained;  however,  in  the  real  world,  there  will  always  be  errors.  As  the  rotation  rate 
increases,  the  error  in  excitation  amplitude  will  increase,  and  depending  on  the  feedback  loop 
coefficients,  the  feedback  loop  may  have  small  error  oscillation. 
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C.  The  CEFP 


In  the  Constant  Excitation  Force  Paradigm  (CEFP),  the  forcing  function  amplitude  is 
constant  irrespective  of  the  gyro’s  rotation.  As  seen  from  Equations  5  and  6,  the  force 
(acceleration)  equations  for  the  gyro  system  are  a  set  of  cross-coupled,  second-order,  linear, 
differential  equations.  This  implies  that  as  the  rotation  rate  changes,  the  X-axis  oscillation 
amplitude  changes.  When  there  is  no  rotation,  the  X-axis  amplitude  is  a  maximum.  As  the 
rotation  increases,  the  X-axis  amplitude  decreases  as  the  Y-axis  oscillation  amplitude  increases. 

This  can  be  thought  of  as  an  iterative  process  that  tracks  the  rotation  rate. 

Thus,  Equation  9  cannot  be  used  to  accurately  determine  the  rotation  rate  because  the 
amplitude  ofx(i)  is  a  function  of  the  rotation  rate  Q.  When  the  gyro  is  rotated  along  the  Z-axis, 
the  solution  to  the  two  equations  can  be  thought  of  as  a  convergent  spiral  of  interactive  forces 
being  generated  between  the  Y-axis  Coriolis  produced  force  and  the  X-axis  Coriolis  produced 
force.  A  forcing  function  produces  a  given  x(t )  which  produces  a  Y-axis  Coriolis  force  which, 
in  turn,  produces  a  X-axis  Coriolis  force  which  changes  the  net  X-axis  force,  changing  x(t) , 
which  then  changes  the  Y-axis  Coriolis  force,  and  so  forth. 

Since  the  electrostatic  forcing  function  is  a  sinusoid,  both  x(t)  and  y(t )  will  be 
sinusoids  at  the  same  frequency,  and  only  the  amplitudes  and  relative  phases  of  x(t)  and  y(t ) 
need  to  be  determined.  Defining  the  X-axis  acceleration  function  as 

Fe  ( t )  =  2 mCly  -  Fx  sin(<a?)  =  FEo  sin(<af ) .  (10) 

It  has  been  shown  in  Equation  4  that  if  the  X-axis  position  is  defined  as 

x  =  x0sin  (cot +  </>),  (11) 


then  a  solution  to  Equation  5  is  when 


x 


=  FEoKa 


(12) 


where  Kco  is  the  physical  system  gain.  Since  co  is  at  the  system  resonance  frequency,  (p  ~  -90 
degrees  and 


x  =  x0a>cos(o)t  -90)  =  -x0rysin(<af). 
First,  define  the  Y-axis  force  as 

-  2 mDx  =  2mQ.caxo  sin(rA) , 


(13) 


(14) 


and  the  Y-axis  position  as 


y  -  y0  sin(ntf  +  9) . 


(15) 
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This  results  in  a  solution  of  Equation  6  where, 


y0=2mQ.coKax0.  (16) 

Again,  since  co  is  at  the  system  resonance  (0  ~  -90  degrees),  the  results  are 

y  -  y0co cos(<stf  -  90)  =  -y0a> sin(ntf) .  (17) 

Revisiting  Equation  5,  the  results  are 

2mCly-Fx  sm(o)t)  =  FEo  %m(o)t)  =  -{imOojf  xr/Kr<)  sin (<x>t)-Fx  %m(o)t)  =  r°  ,  (18) 

Kco 

which  can  be  rewritten  as 


~Fx  =^*-  +  {lm£l(0fxoKa 
or 

-FzK„=x&  +  (2mnaKj) 

or 

-  FyKn 

°  (l  +  (2/nQ  coKC0  )2 ) 

And  thus,  the  equation  for  y0  becomes 

v  ~FxK,(2  mClmK) 

°  [\+{2mO.0)Kj) 

Figure  1  shows  the  nonnalized  trend  for  the  two  signals  x0  and  y„  given  that 
(2mQ.a>Keo)  =  1  at  Q  =  2,000  degrees  per  second. 


(19) 


(20) 


(21) 


(22) 
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Figure  1.  Graph  of  Normalized  Trend  for  Two  Signals 


The  equation  for  x0  predicts  some  interesting  results  when  the  gyro  is  rotated.  For  the 
excitation  or  X-axis: 

•  The  excitation  amplitude  starts  off  at  its  maximum  and  decreased  as  Q  increases. 

•  The  amplitude  decreases  in  a  parabolic  like  fashion  until  it  reaches  Vi,  when 
(2mD.ojKm )  =  1.  In  our  gyro  system,  (2 mQ.coKa)  =  1,  theoretically,  at 
approximately  Q  =  2,000  degrees  per  second. 

•  After  (2mQojKo)  )  =  1,  the  amplitude  starts  decreasing  more  as  the  square  of  Q. 

For  the  sense  or  Y-axis: 

•  The  sense  amplitude  starts  off  at  zero,  and  the  rate  of  increase  decreases  in  a 
parabolic  fashion. 

•  For  small  £2,  the  change  is  almost  linear. 

•  After  (2inQojKio)  =  1,  the  Coriolis  amplitude  decreases  more  linearly  as  £2 

increases.  (The  gyro  was  not  spun  faster  than  2,000  degrees  per  second  to  record 
data  due  to  noise  in  the  communication  channel.) 


The  system  implications  are: 

•  Use  the  Coriolis  signal  to  measure  small  rotation  rates  and  switch  to  the  excitation 
signal  to  measure  larger  rotation  rates. 

•  Use  both  to  measure  rates  in  a  band  in  between. 

•  The  signal  saturation  on  the  sense  output  is  not  caused  by  electronics  saturation. 

IV.  GYRO  EXCITATION  AND  CORIOLIS  MEASUREMENTS 

The  previous  reference  phenomenon  has  been  observed  in  gyro  output  signals.  The 
amplitude  of  the  excitation  signal  starts  decreasing  in  small  steps,  and  as  the  rotation  rate 
increases,  the  size  of  the  steps  increase.  The  sense  signal  starts  increasing  linearly  and  seems  to 
be  reaching  an  asymptote  as  the  rate  increases.  For  example,  in  the  following  test  consisting  of 
uniform  rate  step  size  changes  shown  in  Figure  2,  the  corresponding  sense  signal’s  step  size 
decreases  to  a  point  of  being  difficult  to  notice,  that  is,  the  sense  signal  output  “saturates,”  as 
shown  by  Figure  3.  The  excitation  signal  shown  in  Figure  4  also  decreased  in  distinct  step  as  the 
rate  increases.  At  first,  the  decreases  are  relatively  small,  but  at  larger  rates,  the  decreases  are 
larger  for  the  same  size  rate  increase. 
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Rate  Plot 


Figure  2.  Gyro  Rotation  Rate  in  Degrees  Per  Second 


Figure  3.  Measured  Sense  Signal  Output 
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Figure  4.  Measured  Excitation  Signal  Amplitude 

V.  TEST  DATA  OF  GYRO  CHIP  TO  MEASURE  ACCELERATION 

In  later  tests  on  the  Serial  Number  (SN)  854  gyro,  the  reference  accelerometer  failed,  so 
there  was  no  way  to  determine  the  acceleration  being  applied  to  the  gyro.  Figure  5  shows  a 
portion  of  an  acceleration  test.  The  graph  indicates  the  applied  acceleration  frequency  as  a 
function  of  time,  starting  with  a  15-hertz  signal,  progressing  to  0  hertz,  continuing  with  a  20- 
hertz  signal,  and  so  forth.  Figure  6  shows  the  corresponding  response  on  the  sense  axis 
acceleration  channel.  The  data  sampling  rate  was  such  that  any  acceleration  frequency  greater 
that  approximately  5  hertz  would  be  aliased.  Thus,  all  of  the  acceleration  measurements  in  this 
plot  are  aliased,  and  no  clean  sine  wave  signals  can  be  recorded.  When  the  reference 
accelerometer  worked,  it  measured  low-frequency  acceleration  at  approximately  0.05  g.  The 
measurement  of  the  applied  acceleration  can  be  seen  in  the  changes  in  amplitude  of  the  sense 
axis  channel  corresponding  to  the  time  when  the  acceleration  was  applied. 
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Figure  5.  SN854  Gyro  Expanded  View  Showing  Applied  Acceleration  Frequency 

with  Respect  to  Time 


Figure  6.  SN854  Gyro  Expanded  View  Showing  Acceleration  Response  on  the  Sense 

Acceleration  Channel 
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VI.  CONCLUSION 


To  compensate  for  the  effects  of  vibration  on  the  gyro  required  the  ability  to  measure  the 
vibrations.  Two  methods  to  measure  vibration  were  pursued:  separate  accelerometers  and  the 
gyro  (also  acting  as  accelerometers).  The  gyro  (also  acting  as  accelerometers)  proved  to  be  the 
easiest  to  implement.  Gyro  assemblies  were  vibrated  in  all  three  axes,  and  both  excitation  and 
sense  signal  outputs  were  recorded.  The  gyros  showed  a  white  noise  type  of  reaction  to 
low-frequency  vibration.  At  frequencies  around  co,  standing  waves  were  observed.  True 
environmental  vibrations  would  not  be  at  a  pure  frequency;  therefore,  these  standing  waves 
would  not  occur.  No  feasible  compensation  algorithms  were  developed. 

Equations  were  developed  that  predict  the  changes  in  both  the  excitation  and  sense  signal 
levels  as  rotation  rates  increase  to  higher  rates.  These  effects  are  practically  noticeable  at 
rotation  rates  greater  than  1,000  degrees  per  second.  These  equations  predict  the  saturation  of 
the  sense  output  at  rotation  rates  that  depend  on  the  physical  parameters  of  the  gyro  chip.  For  the 
current  gyro  design,  the  saturation  rotation  rate  was  calculated  to  be  approximately  2,000  degrees 
per  second.  There  has  been  no  testing  above  2,000  degrees  per  second  because  of  the  noise  on 
the  signal  lines  through  the  slip  rings  on  the  high-speed  rate  table  located  at  Aviation  and  Missile 
Research,  Development,  and  Engineering  Center  (AMRDEC).  The  predicted  results  were 
observed  in  rotation  rate  tests.  These  equations  also  suggest  possible  algorithms  to  measure  the 
rotation  rate  while  the  sense  signal  is  approaching  saturation  and  possibly  up  to  twice  the  sense 
axis  saturation  rotation  rate. 
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LIST  OF  ACRONYMS  AND  ABBREVIATIONS 


AMRDEC 

CEAP 

CEFP 

g 

IQ 

MBARS 

mCAS 

MEMS 

SN 

SOI 

co 


Aviation  and  Missile  Research,  Development,  and  Engineering  Center 
Constant  Excitation  Amplitude  Paradigm 
Constant  Excitation  Force  Paradigm 
Gravity 

Inphase  and  Quadrature 
MEMS-Based  Angular  Rate  Sensors 
MicroControlled  Array  Sensors 
Micro-ElectroMechanical  Systems 
Serial  Number 
Silicon-on-Insulator 
Excitation  Frequency 
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APPENDIX 

DERIVATION  OF  THE  EQUATIONS  FOR  THE 
IQ  DEMODULATION  METHOD 


Derivation  of  the  Equations  for  the  In-phase  /  Quadrature  Demodulation 

Method 


This  Appendix  presents  the  demodulation  method  used  to  determine  the  excitation  signal 
amplitude  (at  2rn),  the  excitation  axis  acceleration  (at  co),  and  the  sense  signal  amplitude  (at  to). 
The  method  used  is  referred  to  as  In  phase  /  Quadrature  (IQ)  demodulation.  It  is  easy  to  perform 
in  the  digital  domain.  Starting  with  the  theory,  in  IQ  demodulation  at  co,  the  signal  is  multiplied 
by  both  sin(cot)  and  cos(cot),  summed  over  an  integral  number  of  periods,  squared,  and  added 
together,  and  then  the  square  root  is  taken.  The  following  equations  illustrate  this  process,  where 
N  is  any  positive  integer  and  T  is  the  period  of  co. 
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Along  the  same  line  of  reasoning,  the  equation  for  iz  (from  equation  33  of  the  Excitation 
Appendix)  can  be  written  as 
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Now  define  XIQ  as 
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Now  let  us  assume  that  we  are  computing  IQ  using  a  computer.  Then  the  first 
simplification  is  to  multiply  the  signal  by  a  unity  amplitude  square  wave  at  co  and  2co.  The  sine 
wave  would  be  replaced  with  a  square  wave  that  goes  positive  at  time  0.  The  cosine  wave  would 
be  replace  with  a  square  wave  that  is  delayed  by  90°  or  !4  cycle.  Then  we  notice  that  we  do  not 
have  to  multiply  the  signal  by  the  square  wave  amplitude.  We  only  have  to  sum  the  signal 


A-l 


samples  in  each  time  interval  when  the  square  wave  is  a  constant.  Then  we  would  subtract  the 
integral  sum  when  the  square  wave  amplitude  was  - 1  from  the  integral  sum  when  the  square 
wave  amplitude  was  +1 .  Now  we  notice  that  to  demodulate  at  to,  if  we  divide  the  signal  period 
into  4  equal  time  intervals  (Tl,  T2,  T3,  and  T4)  and  form  four  sums,  i.e.,  the  sum  of  all  the  signal 
measurements  in  each  time  interval,  we  can  calculate  I  and  Q  with  the  same  data.  That  is:  I  =  Tl 
+  T2  -  T3  -  T4,  and  Q  =  -  Tl  +  T2  +  T3  -  T4.  We  now  notice  that  if  we  divide  the  signal  period 
into  eight  equal  time  intervals  and  combine,  i.e.,  add  or  subtract,  the  measurements  in  the  correct 
order,  we  can  demodulate  signals  at  both  co  and  2co  at  the  same  time  using  the  same  signal  sums. 
This  is  the  method  used  to  demodulate  the  excitation  signal  (at  co  and  2co)  and  the  sense  signal  (at 
co).  To  calculate  the  DC  value  or  average  of  a  signal  we  only  need  to  sum  the  signal  in  all  of  the 
time  intervals. 
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